Response to sub-threshold stimulus is enhanced by spatially heterogeneous activity 
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Sub-threshold stimuh cannot initiate excitations in active media, but surprisingly as we show 
in this paper, they can alter the time-evolution of spatially heterogeneous activity by modifying 
the recovery dynamics. This results in significant reduction of waveback velocity which may lead 
to spatial coherence, terminating all activity in the medium including spatiotemporal chaos. We 
analytically derive model-independent conditions for which such behavior can be observed. 

PACS numbers: 05.45.-a,87.18.Hf,05.45.Gg,87.19.Hh 
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Occurrence of spatiotemporal patterns is a generic fea- 
ture of a wide variety of natural systems that are ex- 
citable [l|. Examples of such patterns include travel- 
ing waves in the Belousov-Zhabotinsky reaction 12,] and 
propagating fronts of surface catalytic oxidation [3| . Ex- 
citability in these systems is characterised by a thresh- 
old. A supra-threshold stimulus causes a transition from 
a quiescent to an active state, generating an action po- 
tential (AP). Following this excitation, the system slowly 
returns to rest. During this recovery period the system 
is either fully or partially insensitive to another identical 
stimulus. The dynamical consequences of these proper- 
ties result in distinctive patterns, such as rotating spiral 
waves, that may in turn lead to spatiotemporal chaos ^- 
|6|. While the threshold is a key parameter governing 
excitable systems, the demonstration of stochastic reso- 
nance (SR) [7| and coherence resonance (CR) [8] in such 
media suggest that weak sub-threshold signals could also 
have a significant effect on their dynamics [9, 10]. SR-like 
response resulting from chaotic dynamics in simple sys- 
tems [ll| raises the intriguing possibility that spatially 
heterogeneous activity may enhance the response of an 
excitable medium to sub- threshold signals. 

Understanding how weak signals affect spiral waves 
and other spatiotemporal activity in excitable media is 
especially important because such patterns have criti- 
cal functional consequences for vital biological systems. 
For example, rotating vortices in cardiac tissue that can 
lead to spiral chaos underlie many arrhythmias, i.e., life- 
threatening disturbances in the natural rhythm of the 



heart p^, Il3|. Thus, controlling irregular activity in ex- 



citable media is not only a problem of fundamental in- 
terest in the physics of nonlinear dynamical systems but 
also has potential clinical significance [IJ, ll5|. Exist- 
ing methods of spatiotemporal chaos control in excitable 
systems are almost exclusively dependent on using supra- 
threshold signals, either through a local high-frequency 
source 16|, |l7| or using a spatially extended array [18| . 
Controlling spatial patterns with sub-threshold stimula- 
tion would not only utilize new physical principles, but 
also avoid many of the drawbacks in previously proposed 
schemes. 

In this paper we show that sub-threshold stimulation. 



while having no significant effect on a quiescent medium, 
can induce a remarkable degree of coherence when ap- 
plied on a system with spatially heterogeneous activity. 
Synchronizing the state of activation of all excited re- 
gions ensures that they return to rest almost simultane- 
ously, in the process completely terminating activity in 
the medium. Thus, control of spatially extended chaos 
is achieved efficiently using a very low-amplitude signal. 
We explain the mechanism of this enhanced coherence 
in terms of the role played by sub-threshold stimulus in 
increasing the recovery period of the medium. It signif- 
icantly reduces the propagation velocity of the recovery 
front, thereby increasing the extent of the inexcitable re- 
gion in the medium. We present a semi-analytical deriva- 
tion of the relation between strength and duration of the 
globally applied sub-threshold signal necessary for com- 
plete elimination of spatially heterogeneous activity in 
excitable media. 

A generic model for describing the spatiotemporal dy- 
namics of several biological excitable systems has the 
form: 



dV ^ -hon{V,gi) + l{t) 
dt C„, 



DVV, 



(1) 



where V (mV) is the potential difference across a cellu- 
lar membrane, Cm{= l/iFcm~^) is the transmembrane 
capacitance, D is the diffusion constant ( = O.OOlcm^s"^ 
for the results reported in the paper), /^o^lMAcm" ) is 
the total current density through ion channels on the cel- 
lular membrane, and gi describes the dynamics of gating 
variables of different ion channels. The spatially uniform 
external signal, applied at all points of the simulation do- 
main, is represented by the time-dependent current den- 
sity, I{fiAciR~ ). The specific functional form for lion 
varies for different biological systems. For the results re- 
ported here, we have used the Luo-Rudy I (LRl) model 
that describes the ionic currents in a ventricular cell [19| . 
For all our simulations, the maximum K^ channel con- 
ductance Gk has been increased to 0.705 mS cm~^ to 
reduce the duration of the action potential (APD) [20[. 
To study the effect of sub-threshold stimulus on a sta- 
ble spiral and on spatiotemporal chaos, we have used the 
maximum Ca^+ channel conductance Goj — 0.04 and 



0.05 mS cm~^, respectively. We have explicitly verified 
the model-independence of our results by observing simi- 
lar effects in other realistic channel-based descriptions of 
the ionic current, such as the TNNP model [21|,|22|. 

We consider in turn the response of a single cell, a 1- 
dimensional cable and a 2-dimensional sheet of excitable 
units to a sub- threshold current /. The spatially ex- 
tended systems are discretized on a grid of size L (for 
1-D) and L x L (for 2-D). For most results reported here 
L = 400, although we have used L upto 1200. The space 
step used for all simulations is Sx = 0.0225 cm, while the 
time-step St = 0.05 ms (for 1-D) or 0.01 ms (for 2-D). 
The equations are solved using a forward Euler scheme 
with a standard 3-point (for 1-D) or 5-point (for 2-D) 
stencil for the Laplacian describing the spatial coupling 
between the units. No-flux boundary conditions are im- 
plemented at the edges. The external current is applied 
globally, i.e., I{t) = / in Eq. ([T]) at all points in the sys- 
tem for the duration of stimulation, r. A stimulus {/, r} 
is sub-threshold if it does not generate an action potential 
when applied on a quiescent medium. 

Fig. [1] shows that when sub-threshold stimulation is 
applied to an excitable medium with spatially heteroge- 
neous activity, viz., either a single spiral wave (a-c) or 
spatiotemporal chaos (d-f), there is a striking change on 
the subsequent dynamics of the system. Within a short 
duration (comparable to the APD) there is complete sup- 
pression of all activity in the medium, although in ab- 
sence of this intervention, the existing dynamical state 
would continue to persist for an extremely long time. 
This result is surprising as the weak sub-threshold sig- 
nal appears to be incapable of significantly altering the 
dynamics of an excitable system. To understand this ap- 
parent paradox, we first note that the sub-threshold stim- 
ulation rapidly decreases the number of cells that can be 
excited by existing activity in the medium (Fig. [21 a). 
Indeed, global suppression of activity results when, by 
the end of the stimulation, the number of cells suscep- 
tible to excitation is insufficient to sustain the activity. 
This decrease in their number is because cells tend to 
remain in the recovering state for a longer period in the 
presence of a sub-threshold stimulus. We can see this 
clearly in the response of a single excited cell to a subse- 
quent sub-threshold current / applied for a fixed duration 
r (Fig. [2lb). Increasing / significantly alters the recov- 
ery period resulting in a change in the APD. The time 
f (measured from the initiation of the AP) at which the 
sub-threshold stimulation begins also affects the response 
of the cell to the signal. These results clearly indicate 
that the dominant effect of a sub-threshold stimulus is to 
increase the time period that a cell spends in recovering 
from prior excitation. 

In a spatially extended system, this enhanced recov- 
ery period of the cells results in altering the propagation 
characteristics of the traveling waves. Fig. [3] (a) shows 
a spiral wave propagating in a two-dimensional medium. 



(b) T = 60ms (c) T= 110 ms 




FIG. 1: Pseudocolor plots of transmembrane potential V 
for the two-dimensional LRl model (L = 400) showing the 
elimination of all activity on applying a sub-threshold current. 
The current / is switched on at T = for a duration r = 60 
ms on (a-c) a single spiral, with / =1.6/iA cm~^, and (d- 
f) a spatiotemporally chaotic state, with / =1.8/iA cm~^. 
By T — 150ms, excitation has been effectively terminated 
throughout the simulation domain. 



where each turn of the wave is a region of excited cells, 
with the successive turns separated by recovering regions. 
As the state of the cells evolve with time, it is mani- 
fested in space as movement of excitation and recovery 
fronts. Their propagation speeds are referred to as wave- 
front {cf) and waveback (C5) velocities, respectively. In 
the absence of any external stimulation, Cf — C5, ensur- 
ing that the width of the excited region remains approxi- 
mately constant as the waves travel through the medium. 
However, on applying a sub- threshold stimulus, the wave- 
back velocity becomes significantly lower than that of the 
wavefront which is almost unchanged. Fig. [3] (b) shows 
that, once stimulation begins, C5 quickly decreases to a 
minimum value dependent on /. It then gradually rises 
to eventually become equal to c/ again. For a large sub- 
threshold stimulus /, the waveback velocity rapidly falls 
to its lowest value and changes very slowly thereafter. 
Under these conditions, we can ignore the time- variation 
of C5 for small r and use the time- averaged value C5(/). 
Increasing / leads to an increased difference in the ve- 
locities of the excitation and recovery fronts, c/ — Ch{I) 
(Fig. [31 c). For short stimulus durations, this difference is 
almost independent of r. A significantly lower waveback 
velocity results in the inexcitable region between the ex- 
citation and recovery fronts of a wave becoming extended 
through the course of the stimulation (compare the pro- 
files of APs in a 1-D cable shown in Fig. [4]). This increases 
the overall area of the medium that cannot be excited, 
thereby making it progressively unlikely for the system 
to sustain recurrent activity. 
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FIG. 2: (a) The fraction of cells that can be potentially ex- 
cited (i.e., with V < —60mV) decreases with time during 
stimulation by sub-threshold current / in a two-dimensional 
system with spatiotemporal chaos. The total duration of the 
external signal is r = Ton — Toff = 38 ms. Results shown cor- 
respond to failed (squares) and successful termination (dia- 
monds) of activity in the medium, (b) Effect of sub-threshold 
current on APD of a single cell. The current / is varied keep- 
ing the duration r fixed (= 38 ms). In all cases, the APD is 
shown as a function of the time interval t' between the initia- 
tion of AP and Ton • The inset shows the corresponding effect 
on the AP profile. 



This is explicitly shown for a 1-dimensional cable in 
Fig. |4j When two successive waves propagate along the 
cable, globally applying the sub- threshold stimulus re- 
duces the excitable gap between the recovery front of 
the leading wave (whose velocity C5 has decreased) and 
the excitation front of the following wave (whose veloc- 
ity Cf is unchanged). For a high sub-threshold / applied 
for a long enough duration, the waveback of the first 
wave slows down sufficiently to collide with the succeed- 
ing wavefront. This collision results in termination of the 
excitation front for the second wave which subsequently 
disappears from the medium. 

The above physical picture is fundamentally un- 
changed for a rotating spiral wave with multiple turns as 
shown in Fig. [3] (a). We now use this to propose a sim- 
ple semi-analytical theory for the mechanism by which 
the sub-threshold stimulus suppresses spatially hetero- 
geneous activity. In the absence of any external stim- 
ulus, the width of the excited region of a wave lying 
between its excitation and recovery fronts is / = CfTr^ 
where r^ is the period for which the active cells remain 
excited. This time period is operationally measured as 
the duration for which the transmembrane potential of a 
cell (V) remains above its excitation threshold. On ap- 
plying a sub-threshold external current /, Cf is almost 
unchanged but the resultant waveback velocity, ci){I^t)^ 
which varies with time over the duration of the stimulus 
r, is seen to decrease with increasing /. If / is large or r is 
small, the time- variation of C5 can be neglected and it is 
reasonably well- approximated by the time-independent 
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FIG. 3: (a) Pseudocolor plot of a spiral wave indicating the 
wavelength A, the wavefront velocity c/ and the instantaneous 
waveback velocity Cb(= c/) in the absence of external stimu- 
lation, (b) Time-evolution of Cb during external stimulation 
of duration r = 100 ms using two different / corresponding to 
failed (squares) and successful termination (circles) of activity 
in the medium. During the course of stimulation, c/ is un- 
changed, (c) The average waveback velocity q, (broken lines) 
reduces with / in contrast to c/ (solid lines) 
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FIG. 4: Spatiotemporal evolution of two successive waves 
propagating along a 1-dimensional cable of excitable cells 
when a sub-threshold stimulus is applied over the time in- 
terval [Ton^Toff]. The propagation the wave W2 is blocked 
by the recovery front of Wi due to reduction of its veloc- 
ity, Cfo. The wavefront velocity c/ of both waves are almost 
unchanged. The spatial interval between the two successive 
waves at a time instant is indicated by A. 



average value 05(1) over the stimulus duration. Thus, 
the width of the excited region of the wave increases to 
/(/) = / + [cf — Q)(/)]r. If A is the distance between exci- 
tation fronts of two successive waves in the medium, then 
collision between the recovery front of the leading wave 
and the excitation front of the following wave takes place 
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FIG. 5: Strength-duration curves for a two-dimensional 
medium with the external current applied at all points of a 
quiescent medium (triangles) or a medium with existing ex- 
citation activity, either a single spiral (circles) or spatiotem- 
poral chaos (diamond). The theoretical prediction given in 
Eq. ([2]) is also shown (square). Each (r,/) point is averaged 
over 10 initial conditions. 



when /(/) > A. Thus, for a sub-threshold stimulus /, the 
shortest stimulus duration Tmin necessary to eliminate a 
source of recurrent activity such as a spiral wave is, 



A — CfTr 
Cf -Cb{I)' 



(2) 



Eq. ([2j) provides us with an analytical relation between 
the stimulus magnitude and its minimum duration nec- 
essary for terminating activity in the medium in terms 
of measurable dynamical characteristics of the system. 
Fig.[5]shows that this theoretical strength-duration curve 
for the external stimulation necessary to terminate activ- 
ity matches very well with the empirical data obtained 
from numerical simulations for both single spiral wave 
as well as spatiotemporal chaos. In general, the weaker 
the sub-threshold current, the longer it has to be applied 
in order to alter the dynamical behavior of the system. 
However, there is a lower bound for / below which there 
is no discernible effect of the sub-threshold stimulation 
regardless of its duration. Note that, for values of / just 
above this lower bound, the required Tmin is extremely 
long and the temporal variation of C5 over the duration 
of the stimulation can no longer be neglected. By ex- 
plicitly considering the time-dependence of C5 in Eq. (J2j), 
one can theoretically estimate the value of / where the 
strength-duration curve becomes independent of r. 

The mechanism of the sub-threshold response of ex- 
citable media proposed here depends only on the recov- 
ery dynamics of the system. In detailed ionic models of 
biological excitable cells, this dependence is manifested 
as a decrease in the potassium ion channel conductance 
responsible for the slow, outward K^ current during the 
sub-threshold stimulation. Thus, simplistic models of ex- 
citable media which do not incorporate the effect of exter- 



nal stimulation on the recovery dynamics are inadequate 
to reproduce this enhanced sub-threshold response re- 
ported here. Our results provide a framework for explain- 
ing earlier experimental observations that, in the human 
heart, sub-threshold stimulation not only prevents sub- 
sequent activation 1231 but also terminates certain types 
of arrhythmia J2J, |25| . 

The results reported here may have potential signif- 
icance for controlling spatiotemporal dynamics in sev- 
eral practical situations involving excitable media, such 
as, during clinical treatment of life- threatening arrhyth- 
mias. Current control methods primarily use large supra- 
threshold stimulation to simultaneously activate all ex- 
citable cells. Thus, regions rendered inexcitable through 
prior activity are unaffected, and can subsequently be 
re-activated by any remaining excitation after the exter- 
nal stimulation is removed, leading to failure of control. 
By contrast, the sub-threshold stimulation method de- 
scribed here slows the recovery of excited cells, thereby 
reducing the pool of cells available for excitation by ex- 
isting activity in the medium. As this approach, in gen- 
eral, requires lower energy, it suggests a complementary 
approach for efficiently terminating spatially extended 
chaos in excitable systems. The proposed mechanism 
may also provide a key to understand how spatially ir- 
regular activity in biological systems (e.g., fibrillation) 
can be controlled by signals strongly attenuated during 
passage through the intervening medium [26|, l27| . 
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